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Abstract 


The  requirement  for  a  highly  efficient  plasma  source  that  can  operate  in  an  electric 
propulsion  capacity  has  led  to  renewed  research  on  helicons.  Helicons  have  long  been 
known  to  be  a  highly  efficient  (>90%)  and  high  density  (>1013  cm'3)  method  to  generate 
plasma.  It  has  not  been  until  recently  though,  that  these  sources  have  been  given  interest  by 
the  propulsion  community.  Much  recent  research  toward  this  end  has  focused  on 
detachment  of  ions  from  the  magnetic  field  lines  which  can  be  as  great  as  a  few  hundred 
Gauss.  Investigations  to  date  have  shown  that  this  may  be  possible  through  ambipolar 
effects  or  acceleration  through  a  double  layer.  Preliminary  investigations  suggest  electrons 
that  are  generated  near  the  diverging  section  of  the  axial  magnetic  field  will  enhance  the 
ambipolar  acceleration  of  ions  out  of  the  device.  Toward  this  end,  a  radio  frequency  (RF) 
antenna  that  propagates  a  wave  toward  the  converging  section  of  the  device  (diverging 
section  of  the  magnetic  field)  will  dictate  where  the  bulk  of  the  RF  energy  is  absorbed.  This 
has  led  to  the  design  of  a  conical  antenna.  Preliminary  work  based  on  a  uniform  density 
model  suggest  that  at  smaller  radii  there  is  an  increase  in  electron  density  and  that  there  are 
further  ways  to  control  the  density  and  temperature  based  on  considerations  such  as 
pressure,  frequency,  and  antenna  length.  The  focus  of  this  investigation  will  be  to  observe 
and  compare  the  differences  in  the  wave  fields  (Br,  B0,  Bz)  between  a  traditional  cylindrical 
helicon  and  that  of  a  conical  helicon. 
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Nomenclature 


a 

Bo 

Br 

Be 

B-_ 

E 

e 

f 

g 

i 

Jm 

j 

k 

k 

m 

n 

n0 

ne 

p,r 

r 

t 

Y,,m 

x,y,z 


fixed  radius 

static  magnetic  field 

radial  component  of  magnetic  field 

azimuthal  component  of  magnetic  field 

axial  component  of  magnetic  field 

electric  field  wc  =  cyclotron  frequency 

electric  charge  cop=  plasma  frequency 

driving  frequency 

generic  function 

imaginary  number 

Bessel  function 

current  density 

wave  number 

wave  vector 

mode  number 

density  (plasma  +  neutral) 

plasma  density 

electron  density 

associated  Legendre  polynomial 

radius 

time 

spherical  harmonic 
Cartesian  units 


Ho  =  permeability  of  free  space 
(f>  =  spherical  angle 
6  =  cylindrical/spherical  angle 
Gwan  =  fixed  angle  of  cone 
co  =  driving  frequency 


I.  Introduction 

The  term  ‘helicon’  was  first  used  in  1960  by  Aigrain  [1],  to  describe  waves  that  propagate  in  the  presence  of  a 
magnetic  field  for  the  frequencies  between  the  ion  and  electron  cyclotron  frequencies.  Helicon  waves  are  a 
subclass  of  whistler  waves  which  were  first  investigated  by  Storey  in  the  1950’s  [2],  Whistler  waves  are  right- 
handed  circularly  polarized  waves  that  propagate  in  free  space  differing  from  helicons  in  two  main  regards:  (1) 
helicons  are  of  such  a  lower  frequency  than  whistler  waves  that  the  electron  gyrations  (being  so  much  faster)  may  be 
neglected  and  (2)  helicons  propagate  in  a  bounded  system  [3],  specific  to  this  study,  a  cylindrical  or  conical 
boundary.  Helicon  waves  were  initially  studied  in  metals  and  semiconductors  by  Libchaber  [4]  and  Rose  [5]  in  the 
early  1960’s  who  observed  the  wave  propagation  at  frequencies  much  less  than  the  electron  cyclotron  frequency  and 
whose  wave  vector  k  was  parallel  to  the  magnetic  field.  Soon  afterwards,  helicons  were  observed  in  gaseous  plasma 
by  Lehane  and  Thoneman  [6],  Kuckes  [7],  and  Nyack  and  Christiansen  [8]  who  all  observed  right-handed  circularly 
polarized  waves  in  a  laboratory  experiment.  The  experiments  were  run  from  power  levels  of  several  hundred  watts 
to  three  kilo-watts  while  a  dc  solenoidal  field  of  several  hundred  Gauss  was  applied.  All  three  found  comparisons  to 
agree  well  with  the  theory  of  plasma  wave  propagation  in  an  insulated  cylinder  for  a  uniform  plasma  profile 
originally  derived  by  Klozenberg,  McNamara,  and  Thoneman  [9]  (KMT  theory).  However,  the  experiments  were 
limited  in  scope  by  the  diagnostics  and  equipment  at  hand  and  it  wasn’t  until  the  1970’s  when  Boswell  presented  his 
work  on  waves  in  gaseous  plasma  [10-13]  and  Chen  [3,  14-27]  began  to  investigate  helicons  that  much  experimental 
data  began  to  be  produced.  Thus,  the  great  majority  of  theory  and  experimental  data  that  is  relevant  to  this  study  was 
done  by  Boswell,  Charles,  and  Chen.  Specifically,  the  resonant  absorption  of  waves  in  plasma  at  or  near  the  lower 
hybrid  frequency  as  well  as  downstream  density  peaks  and  most  recently  double  layers  is  what  this  work  will 
ultimately  investigate. 
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II.  Cylindrical  Wave  Fields 


The  description  of  helicon  waves  follows  the  analysis  of  Chen  and  others  where  the  basic  equation  to  be  solved  is 
the  simplified  Ohm’s  law 

neE  =  jxB0  (1.1) 


where  the  static  applied  magnetic  field  is  B(j  =  BQZ  and  we  use  Maxwell’s  equations 


VxE 


8B_ 

dt 


VxB  =  /j0j 


VUB  =  0 


vq/  =  o 


with  perturbations  of  the  form  exp  i(mO  +  kz  —  COt)  to  obtain 


B  = 


f  V1 

CO  jUoeno 
k  B0  j 


VxB 


where  we  let 


2 

CO  fJ()enQ  _  °0  60 P 

k  B0  k  cocc 1 


And  equation  ( 1 .6)  then  becomes 


VxB  =  aB 


and  finally  taking  the  curl  of  equation  (1.8)  yields 


V2B  +  a2B  =  0 

Writing  this  in  cylindrical  coordinates  for  the  z-component  gives 


K+-K  + 

r 


v  ~  r2  J 


B=  0 


(1.2) 

(1.3) 

(1.4) 

(1.5) 

(1.6) 


(1.7) 


(1.8) 


(1.9) 


(1.10) 


where  T  =  a~  +  k 
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r  =  5  cm 


r  =  10  cm 


r  =  1 5  cm 


Figure  1.  Evolution  of  Cylindrical  Helicon  Wave  Fields  with  Increasing  Applied  Magnetic  Field. 


There  must  be  a  finite  solution  at  r=0  ->  Jm(Tr),  so  that 

Bz=C,Jm(Tr) 


(1.11) 


The  r  and  #  components  are  then  obtained  from  equation  (1.8) 


Iff! 

—  B  -ikBa  =  aB„ 


(1.12) 


ikB„  -  B  =  aBa 


(1.13) 


These  can  then  be  algebraically  manipulated  into  solving  for  the  r  and  <9  components  of  the  magnetic  field  in  terms 
of  Bessel  functions 


B,  = 


iC\ 


-aJ-kJ. 


\ 


m  m 


(1.14) 
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(1.15) 


Be 


kJm+aJm 


) 


We  can  now  analyze  what  the  field  structure  should  look  like  given  some  basic  input  parameters.  For  the  cases 
considered  here,  we  assume  a  uniform  electron  density  of  1012  cm"3,  13.56  MFlz  driving  frequency,  and  an  m=+l 
helical  antenna.  We  take  a  look  at  three  separate  cases  of  varying  radial  size,  namely,  cylinders  of  5,  10,  and  15  cm 
radii.  The  results  are  indicated  in  Figure  1.  Flere  we  find  that  for  the  density  selected  helicon  waves  will  not 
propagate  at  lower  magnetic  fields  if  the  radius  is  too  large.  At  100  Gauss  external  field,  only  the  5  cm  radius 
geometry  will  support  helicon  waves.  It  is  not  until  the  applied  field  reaches  250  Gauss  that  the  10  cm  geometry 
will  support  wave  propagation  and  until  about  500  Gauss  that  the  15  cm  geometry  will  support  helicon  waves. 
Additionally,  once  the  applied  external  magnetic  field  is  great  enough  to  allow  helicon  wave  propagation,  further 
increase  in  applied  field  does  not  change  the  structure.  At  least  not  until  more  input  power  is  available.  A  higher 
level  of  input  power  will  raise  the  electron  density  and  the  helicon  wave  field  will  start  to  dissipate  until  a  higher 
magnetic  field  is  applied.  Clearly,  there  is  a  delicate  balance  between  the  input  power  and  applied  magnetic  field 
necessary  to  sustain  helicon  wave  propagation. 

Being  able  to  address  the  profile  of  the  radial  magnetic  field  is  what  differentiates  a  purely  inductive  discharge  from 
a  helicon  discharge.  So  if  we  are  to  investigate  different  geometries  such  as  the  conical  one  proposed  here,  we  will 
need  to  perform  a  similar  analysis  to  that  done  in  cylindrical  geometry  to  look  at  the  new  radial  profiles. 


III.  Spherical  Wave  Expansion  for  Conical  Geometry 

In  cylindrical  coordinates,  helicon  waves  are  described  by  perturbations  of  the  form, 

exp  i{m6  +  kz-  cot )  (2. 1) 

where  the  coordinate  system  used  is  r  +  0  +  Z  .  This  represents  waves  traveling  in  the  +z  direction  in  time  and 
rotating  either  clockwise  (m>0)  or  counterclockwise  (m<0)  with  respect  to  6.  Flowever,  when  examining  helicon 
waves  in  a  spherical  coordinate  system  to  describe  wave  propagation  in  a  cone,  the  perturbations  must  take  a 
different  form.  To  begin  the  analysis  we  assume  spherical  perturbations  are  of  the  form 

exp  i(m<j>  +  kr-(ot)  (2.2) 

where  the  coordinate  system  is  r  +  9  +  (f) .  This  represents  waves  traveling  in  the  r  direction  in  time  and  rotating 
either  clockwise  (m>0)  or  counterclockwise  (m<0)  with  respect  to  (/>.  The  perturbations  described  here  assume  wave 
propagation  in  the  r  direction  where  the  cone  wall  is  given  by  a  constant  6  =  0waU  from  the  x-axis.  The  situation  is 
illustrated  below  in  Figure  2. 
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When  utilizing  a  cylindrical  coordinate  system  to  study  helicon  waves,  the  system  inherently  allows  for  the 
simplicity  of  applying  a  constant  magnetic  field  Bn  in  the  z-direction;  or  the  direction  of  wave  propagation. 
However,  when  analyzing  wave  propagation  in  spherical  geometry,  a  constant  magnetic  field  in  the  Cartesian  or 

cylindrical  z-direction  now  must  be  given  in  terms  of  the  different  unit  vectors;  r,  6 ,  (f)  .  The  system  described 
above  is  shown  in  Figure  3  but  now  includes  the  applied  magnetic  field  in  the  Cartesian/cylindrical  +z-direction. 
The  objective  is  to  describe  the  magnetic  field  in  terms  of  the  spherical  unit  vectors. 
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To  describe  the  unit  vectors  in  the  z-direction,  we  must  first  write  the  spherical  unit  vectors  in  terms  of  their 
cartesian  counterparts, 

r  =  sin  6  cos  (fix  +  sin  6  sin  <f>y  +  cos  6z 
6  =  cos  6  cos  (j)x  +  cos  6  sin  tf>y  -  sin  dz  (2.3) 

^  =  -  sin  <f>x  +  cos  <j)y 

These  are  the  normalized  spherical  unit  vectors.  We  can  now  algebraically  manipulate  the  above  3  equations  and 
solve  for  the  unit  vector  in  the  z-direction  in  terms  of  l',0,(f)  which  yields 

z  =  cos  Or  —  sin  66  (2.4) 

Therefore,  when  describing  the  static  magnetic  field  applied  in  the  z-direction,  we  can  relate  the  different  coordinate 
systems  by, 

B  =  B0  z  =  B0  (cos  Or  -  sin  60)  (2.5) 

We  can  now  begin  to  describe  the  wave  propagation  in  spherical  coordinates  beginning  with  Ohm’s  law  as  in 
equation  (1.1), 


neE  =  jxB0  (2.6) 

where  B0  is  given  above  in  equation  (2.5).  Again  we  make  use  of  Maxwell’s  equations 


—  8B 

VxE  = - 

(2.7) 

dt 

VxB  =  /uQj 

(2.8) 

VDS  =  0 

(2.9) 

< 

ll 

o 

(2.10) 

and  perturbations  given  by  exp  i{m(j)  +  kr  —  cot)  so  that  equation  (2.7)  becomes 

VxE  =  icoB  (2.11) 


Taking  the  curl  of  equation  (2.6)  to  obtain 

YxE  VXJ'XB"  (^oC^)-/~(-/W)^o  +  -/(VDgo)--go(VD-/)  (212) 

me  ne 

The  third  and  fourth  terms  of  equation  (2.12)  are  zero  as  a  result  of  Maxwell’s  equations  (2.9)  and(2.10). 
Additionally  the  second  term  is  zero  because  the  applied  magnetic  field  is  constant  in  all  directions.  If  we  now 
equate  the  left  hand  sides  of  equation  (2.1 1)  and  (2.12)  we  arrive  at 
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(2.13) 


icoB  = 


(Wy 


ne 


In  terms  of  the  wave  perturbations,  j  is  similarly  defined  as 

j  =  j  exp  i(nnf>  +  kr  -  cot) 

If  we  now  expand  equation(2.13),  we  obtain. 


(2.14) 


icoB  =  —  - 
ne 


[  50  cos  #r  -  50  sin  ##  |[|  —  r  +  —  -^—0  + 


(b 

dr  r  d6  r  sin  6  d(j)  ) 


v 

1-1 

A 

). 

J  J 

which  upon  simplifying  the  inner  product  becomes, 


r 


icoB  =  — < 
ne 


d 


i  a 


Bn  cos 0 - Bn  sin# - 0 


dr 


r  d6 


Inserting  equation  (2.14)  for  j  into  equation  (2.16) 


icoB  =  —< 
ne 


r 


d 


i  a 


a 


B '  cos# - Bn  sin# - # 


dr 


r  d6 


j  exp  i(nnf)  +  kr  -  cot) 


J 


(2.15) 


(2.16) 


(2.17) 


However,  the  partial  derivative  with  respect  to  0  is  zero  because  the  current  term  is  assumed  to  be  constant  across 
the  radial  arc  from  -0wan  to  +0waii. 

Note:  This  assumption  does  not  include  a  perturbation  in  the  0-direction.  This  assumption  was  based  on  the  helicon 
wave  analysis  done  in  cylindrical  geometry  where  the  basic  assumption  of  a  uniform  density  gives  rise  to  plane 
wave  propagation.  In  spherical  coordinates,  similar  waves  are  assumed  to  propagate  although  the  waves  are 
spherical  waves  of  a  constant  radius  ‘r’. 

Simplifying  equation  (2.17)  will  yield, 


icoB  = 


B0ik  cos  #  ( j  exp  i ( m</>  +  kr  -  cot ) ) 


ne 


(2.18) 


which  becomes 


COS 


0)1 


cone 


Inserting  equation  (2.8)  fory 


B  = 


[B0k  cos#) 


VxB 


jU0cone 


(2.19) 


(2.20) 
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If  we  define  the  constants  in  equation  (2.20)  as 


ju()cone 

B0k 


(2.21) 


Then  equation  (2.20)  becomes 


aB  =  (cos6,)Vx5 


(2.22) 


Taking  the  curl  of  both  sides  of  equation  (2.22) 

aVxB  =  Vx  |(cos6*)  VxZ?j  (2.23) 

The  analysis  becomes  further  complicated  when  considering  the  curl  in  equation(2.23),  there  is  the  additional  cos 0 
term.  Equation  (2.23)  then  becomes 


(2.24) 


And  in  order  to  ultimately  arrive  at  the  second  order  differential  equation  for  the  magnetic  field  structure,  the 
Laplacian  in  spherical  coordinates  is  needed 


v!=il 

r2  dr 


V 


dr 


+  - 


1 


r  sin  6  dO 


3 

do 


+- 


i 


r~  sin"  6  dtp2 


(2.25) 


For  reference,  the  curl  in  spherical  coordinates  is  given  by 


If  we  attempt  to  look  at  this  from  a  different  perspective  and  assume  that  because  the  conical  angle  we  are  looking  at 
in  the  lab  is  small  (~  4° ),  we  can  approximate  B  as 
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(2.27) 


B  =  B0z  =  B0r 


We  will  still  assume  perturbations  of  the  form  expi(m<f>  +  kr  —  (Ot)  and  will  reconsider  the  problem  beginning 
with  equation  (2.13)  and  expanding  such  that 


icoB  =  —  < 
ne 


(v) 


/  8  -  1  d  2  1  3 

0  — r  + - 0  + - 

dr  r  d6  rsin#  d(j) 


<t> 


y 

1-1 

A 

J. 

J  J 

(2.28) 


Substituting  equation  (2.14)  into  equation  (2.28)  yields 


icoB  =  —  I  B()  —  (  /  exp  i(m</>  +  kr  -  cot))  >  = 


ne 


dr 


BJk 


ne 


If  we  now  substitute  Maxwell’s  equation  (2.8)  for  j  into  equation  (2.29)  we  arrive  at 


(2.29) 


B=  —  V  x  #  =  —V x  B 
ju0cone  a 


(2.30) 


where  a  is  defined  identically  as  in  equation  (2.21).  Taking  the  curl  of  both  sides  of  equation  (2.30)  where 
V  x  x  B j  =  V  B  where  the  second  term  is  zero  according  to  Maxwell’s  equation  (2.9).  Equation 

(2.30)  then  becomes 


-V2B  =  aVxB  =  a2  B 
V2B  +  cc2B  =  0 

Which  is  the  Helmholtz  equation  for  spherical  coordinates  and  is  written 

i  u  f  .  ndH)  1  ?2" 

dr  r~  sin0  d0(  dd 


_L  d_ 

r2  dr 


r  2dB ^ 

r 


V 


J 


r  sin#  dQ 


1  d2B 

r2  sin2  #  d(jf 


+  a~B  =  0 


This  equation  is  separable  using  the  approach 

B(r,0,t)  =  R(r)®(0)®(4) 

When  substituting  equation  (2.34)  into  equation  (2.33)  yields 


j _d_ 
R  dr 


V 


dR 

dr 


2  2  1 

+  CTr~  +- 


0sin#  d6 


sin# 


V 


#0 

do 


1  d2®  n 

+ - =0 


Osin-  #  dfr 


(2.31) 

(2.32) 


(2.33) 


(2.34) 


(2.35) 


We  can  take  the  terms  that  only  involve  <|>  of  the  last  term  in  equation  (2.35)  and  set  it  equal  to  a  constant  -m2  we 
obtain 
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This  transforms  equation  (2.35)  such  that  the  last  two  terms  are  now  a  function  of  6  only  and  the  constant  m.  We 
can  set  these  terms  equal  to  another  constant  -1(1+1)  which  will  be  helpful  when  solving  the  radial  components.  The 
#  components  of  equation  (2.35)  then  become 


1 


d 


0sin#  d6 


sin# 


V 


de 

de 


+-^  =  -/(/+ i) 

sin2  6  X  ’ 


(2.37) 


Equation  (2.36)  has  the  familiar  solution  as  a  simple  harmonic  oscillator  while  equation  (2.37)  is  recognized  as  the 
associated  Legendre  equation.  The  solution  to  both  are  given  by 

$>(</))  =  eim<l  (2.38) 

0(#)  =  if  (cos#)  (2.39) 

In  which  1=0, 1,2,...  and  m  runs  over  the  integer  values  from  -1  to  1,  and  if  ( COS # )  are  the  associated  Legendre 

polynomials.  If  1  is  not  an  integer,  the  solution  to  equation  (2.37)  diverges  for  cos9  =  1  or  -1  (9  =  0  or  n).  Often  the 
solutions  of  0  and  <D  are  combined  such  that 


0(#)<D(y)  =  if  (cos#)ew  =  If  (#y) 


(2.40) 


where 

yields 


are  known  as  spherical  harmonics. 


Now  examining  the  radial  components  of  equation  (2.35) 


d_ 

dr 


dR 

dr 


2  2 

a  r 


-/(/  +  l)]i?  =  0 


(2.41) 


Equivalently,  this  can  be  written  as 

r2^-  +  2r  —  +  \a2r2-l(l  +  l)\R  =  0 
dr2  dr  L  V 


Making  the  appropriate  substitution  so  that 


R(r) 


£0 

(. ar )/2 


Transforms  equation  (2.42)  to 


r 


2 


d2Z  dR 

— y-  +  r - + 

dr  dr 


2  2 

a  r 


Z  =  0 
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The  solution  to  this  equation  is  given  by  the  Bessel’s  function  of  1+1/2  order  which  are 


.  Combining  these  with  the  original  substitution  of  equation  (2.43)  yields 


J,+y2(ar ) 

( ar )^2 


Ji+y\ar) and 


(2.45) 


which  are  the  spherical  Bessel  functions  of  the  first  kind.  The  terms  have  been  dropped  to  ensure  that 

the  solution  is  finite  at  the  origin.  If  we  now  recombine  solutions  to  the  spherical  Helmholtz  equation  as  defined  in 
equation  (2.34)  we  obtain 


00  / 

B(r,0,<t>)  =  YJY,  aimJi{ar)Yr  (W)  (2-46) 

/=0  m=—l 

where  «/,„  are  determined  by  the  boundary  conditions. 

Note:  This  is  of  course  the  simplified  version  of  the  traveling  spherical  waves. 

Key  assumptions  made: 

1)  Perturbation  is  of  the  form  exp  i{m(j)  +  kr  —  (Ot)  .  We  know  that  this  is  incorrect  because  it  violates 
conservation  of  energy.  The  energy  density  stored  in  the  wave  is  equal  to  the  square  of  its  amplitude. 
However,  the  outgoing  wave  increases  its  area  by  a  factor  of  r  .  This  would  imply  that  the  energy  of 
the  wavefront  increases  as  r2.  This  cannot  be  true  and  we  can  eliminate  this  fact  by  dividing  the 
amplitude  by  r  (which  divides  the  energy  density  by  r2).  So  the  new  (correct  form)  of  perturbations  to 
the  spherically  traveling  waves  must  be 

exp  iimd)  +  kr  —  cot) 

g^go  - -  (2-47) 

r 

2)  The  magnetic  field  is  given  by  B  =  B()Z  =  B0 r .  As  shown  earlier,  the  actual  applied  field  needs  to 
be  described  in  terms  of  the  spherical  coordinate  unit  vectors  which  yielded, 

B  =B0z  =  Bo(cos0r-sin&d)  (2.48) 

Therefore,  the  next  part  of  the  analysis  will  be  to  correct  for  the  first  assumption. 


This  time  we  utilize  the  expansion  of  equation  (2.13)  for  B  =  B()Z  =  BQr  and  j  given  by 


j  =  Jo 


Qxpi(m<f>  +  kr- cot) 
r 


(2.49) 


Upon  substitution  yields 
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-  1  d 

iojB  =  <  B(  — 

ne  or 


J 


exp  i{ni(f)  +  hr  -  cot) 


Evaluating  the  derivative  product  gives 


ne  v  r 


''  j  exp  i(m(f)  +  kr  -  cot )  =  —  — 

ne  1  r 


r2  J 


r2  J 


Substituting  Maxwell’s  equation  (2.8)  for  the  current  density  j  gives 


B  = 


Bn 


( ik  1  ^ 


i  Hifione  V  / 

Taking  the  curl  of  both  sides  of  equation  (2.52)  yields 

Bn 


VxB 


j 


f 


VxB  = 


A 


yiju0cone  J 


Vx- 


^  ik 

V  r 


r2  J 


VxB 


(2.50) 


(2.51) 


(2.52) 


(2.53) 


At  this  point,  an  analytical  solution  appears  unlikely  and  the  equation  that  will  have  to  be  numerically  solved  is 


ex' k  ' r  B  +  ( ikr  - 1)  V : 


ikr  - 1 


Vx5  =  0 


(2.54) 


J 


The  attempt  of  this  analysis  is  similar  to  that  in  cylindrical  coordinates.  Understanding  the  plasma  magnetic  field 
structure  showed  when  helicon  waves  could  be  supported.  Therefore,  the  same  approach  was  taking  for  analyzing  a 
conical  antenna,  although  ultimately  the  analysis  became  much  more  complex.  At  this  juncture,  the  next  step  is  to 
either  numerical  solve  the  equations  involving  spherical  waves  or  to  assume  plane  wave  propagation  in  spherical 
coordinates  subject  to  the  conical  wall  boundary  conditions  which  greatly  simplifies  the  solution  to  the  spherical 
harmonics  obtained  above. 


IV.  Hardware  and  Initial  Setup 

Current  hardware  consists  of  3  quartz  cylinders,  each  ~  45  cm  length.  Two  are  straight  cylinders  with  3.81  cm  and 
6.35  cm  diameters,  while  the  last  is  a  conical  design  starting  at  3.81  cm  diameter  and  ending  at  10.16  cm  diameter  as 
shown  schematically  below  in  Figure  4. 
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▲ 


10.16cm 


Figure  4.  Schematic  of  conical  quartz  tube. 


Each  quartz  tube  is  mounted  through  a  conflat  flange  that  is  connected  to  a  0.5  m  diameter,  1  m  long  cylindrical 
vacuum  chamber  that  is  pumped  to  a  base  pressure  of  lxlO'6  torr  through  a  250  1/s  turbomolecular  pump  which  is 
sufficient  for  helicon  operation  in  the  10"4-10"'  torr  range.  The  pressure  is  monitored  through  a  convectron  gauge  in 
the  high  pressure  range  (>10  mTorr)  and  by  a  cold  cathode  gauge  in  the  low  pressure  range  (<  10"4  torr).  However, 
in  the  intermediary  range  of  10"4-10"2  where  neither  gauge  gives  the  most  accurate  pressure  a  capacitance  manometer 
calibrated  to  this  pressure  range  is  used.  Both  vacuum  chamber  and  quartz  tube  are  mounted  on  a  vibrationally 
insulated  optics  table. 


Electrostatic  Ion 
Energy  Analyzer 


90  GHz  Microwave 
Interferometer 


0-1200  Gauss 
DC  Magnet  Coils 

/■ 


RF  Compensated 
Langmuir  Probe 


K  7  MM  MM  ^  Hcjlf-Turn 
A _ P  Helical  Antenna 


Quartz  Tube 


iTurbo  &  Roughing 
Pumps 

IZ 


Matching  Network 


- / 

RF  EMI  Shielding 


7 


Oscilloscope 

Power  Meter 

0-100  MHz 
Waveform  Generator 


0.3-35  MHz  0-1000W 
Linear  Amplifier 


Figure  5.  Schematic  of  AFRL  Experimental  Helicon  research  station. 
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Power  to  the  antenna  /  plasma  is  supplied  by  a  0.3-35  MHz  lkW  broadband  RF  generator  fed  through  a  autotune 
matching  network  to  the  m=+l  helical  antenna  made  of  copper  strap.  The  static  magnetic  field  is  supplied  by  4 
water-cooled  magnets  that  can  supply  a  continuous  uniform  field  up  to  1200  Gauss.  A  schematic  of  the  current 
setup  with  chamber,  magnets,  antenna,  matching  network,  and  quartz  tube  are  shown  in 

Current  diagnostics  consist  of  an  RF  compensated  Langmuir  probe  with  a  copper  compensation  electrode  and  4 
chokes  near  the  probe  tip  to  minimize  RF  pickup  at  the  driving  frequency  of  13.56  MHz  and  its  first  harmonic  27.12 
MHz.  The  probe  design  was  modeled  after  that  of  Sudit  and  Chen  [28].  Similar  to  previous  helicon  work,  a  90  GHz 
microwave  interferometer  has  been  used.  Both  diagnostics  are  currently  ready  for  use  upon  proper  RF  shielding  of 
interference. 

Lastly,  photos  of  the  discharge  have  been  taken  to  demonstrate  the  capability  of  generating  RF  and  helicon 
discharges.  Photos  of  inductively  coupled  plasma  are  shown  below  in  Figure  6  for  no  applied  external  magnetic 
field. 


Power  <  100  Watts 


Power  ~250  Watts 


Figure  6.  Plasma  generated  with  m=+l  antenna  with  no  applied  external  B-field. 

Photos  for  what  is  believed  to  be  a  helicon  discharge  at  two  different  power  levels  and  800  gauss  of  external  field 
are  shown  in  Figure  7. 


Figure  7.  Plasma  generated  with  m=+l  antenna  with  800  Gauss  applied  field. 
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V.  Conclusions  and  Future  Work 


This  study  has  revealed  the  fact  that  cylindrical  geometry  inherently  simplifies  the  analysis  toward  arriving  at  the 
helicon  dispersion  relation  and  evaluating  the  wave  structure  as  opposed  to  a  different  geometry  such  as  a  conical 
antenna  analyzed  with  spherical  coordinates.  This  is  evident  by  the  planar  wave  propagation  in  cylindrical  geometry 
and  the  analysis  is  complicated  by  spherical  wave  propagation  and  the  requirement  of  conservation  of  energy. 
However,  the  analysis  has  led  to  a  differential  equation  that  should  be  straightforward  to  analyze  numerically  and 
ultimately  this  is  the  next  step.  The  goal  of  solving  in  spherical  coordinates  is  to  determine  the  wave  structure  and 
then  compare  the  solution  to  that  of  the  cylindrical  helicon.  The  numerical  and  analytical  aspect  of  this  will  be 
complemented  with  experimental  measurements  using  b-dot  probes  to  look  at  the  radial  fields. 

The  projected  work  will  include  comparisons  between  the  three  quartz  tubes  currently  available,  two  cylinders  of 
different  radii  and  the  conical  tube.  To  date,  the  author  is  not  aware  that  anyone  has  attempted  a  conical 
configuration  for  the  antenna  and  insulating  wall  for  plasma  wave  propagation.  Based  on  a  separate  uniform  density 
model,  we  arrived  at  the  simple  expectation  that  a  smaller  radius  tube  will  yield  higher  density  plasma  for  the  same 
power,  magnetic  field,  pressure,  and  frequency.  Consequently,  comparisons  between  the  three  quartz  tubes  as  far  as 
axial  and  radial  density  and  temperature  profiles  will  be  done  through  the  use  of  RF  compensated  Langmuir  probe 
analysis  verified  against  microwave  interferometry  data.  The  above  measurements  will  have  to  be  taken  for 
different  power  levels  and  magnetic  field  configurations  to  determine  the  initiation  of  the  helicon  wave  discharge  as 
described  in  Section  II. 

Later  on  in  this  study,  attention  will  be  given  to  the  observation  of  a  double  layer  by  Charles  at  fields  of  250  gauss 
[29-3 1],  as  well  as  the  effect  of  pressure  on  plasma  potential  from  the  antenna  downstream.  The  majority  of  current 
helicon  plasma  operate  at  pressures  >  1  mTorr  where  the  mean  free  path  for  collisions  is  much  less  than  the  antenna 
and  radii  dimensions.  However,  as  pointed  out  by  Charles,  at  lower  pressures  (0.2-0. 4  mTorr)  the  mean  free  path 
can  be  increased  to  the  order  of  device  dimensions  and  this  is  believed  to  be  what  allows  the  formation  of  the  double 
layer.  Therefore,  the  importance  of  investigating  the  pressure  effects  to  determine  whether  this  is  the  case  should  be 
explored. 
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Helicons:  The  Why? 


•  Well. .  .everyone  else  is  doing  it. . . 


Company /Institution 

Research  Time 

Reference 

AFRL 

2006-Present 

Personal  Involvement 

University  of  Illinois 

2005-Present 

Personal  Involvement[28] 

Starfire  Industries 

2005-Present 

NASA  (Huntsville) 

2006-Present 

[29] 

University  of  Washington 

2004-Present 

[30] 

University  of  Texas 

1998-Present 

[31-36] 

UCLA 

1980’s-Present 

[3,  14-27] 

University  of  Wisconsin 

2003-Present 

[37,  38] 

Stanford  University 

2005-Present 

West  Virginia  University 

2000-Present 

[39,  40] 

Cornell  University 

1991-??? 

[41] 

Europe  and  Japan 

2003?-Present 

[42-44] 

Australian  National  University 

1965-Present 

[10-13,45-53] 

Actually 

•  Electrode-less 

•  High  Density  at  Low  Power 

1019nr3  @  100’s  Watts 

•  High  Ionization  Efficiency  >  90% 

•  Variable  ISP  /  Exhaust  Velocity 

•  Wide  range  of  gas  /  molecular 
propellant 
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What  is  a  Helicon  and  Why? 


Subclass  of  Whistler  Waves. .  .what  is  a  whistler  wave? 

-  Wave  travels  along  magnetic  field  lines  in  free  space 

-  Right  Hand  Polarization 

•  The  electric  field  rotates  clockwise  when  looking  in  the  direction  of  propagation 


long 

whistlers 


Helicon  Wave 

-  Bounded  whistler  wave 

-  Frequency  between  ion  and  electron  gyro-frequency 


Helicon  Discharge 

High  Density,  Axially  Magnetized  Plasma  Column 

cociU  coLH  □  con  COce 
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CCP,  ICP,  Helicon 


CCP  -  Capacitively  Coupled  Plasma 

-  Electrons  respond  to  E-Field  ->  Gain  Energy 

-  Can  ionize  gas  directly  or  indirectly 

ICP  -  Inductively  Coupled  Plasma 

-  Time  varying  E-field  generates  a  time  varying  magnetic  field  in  the 
longitudinal  direction 

-  Induces  azimuthal  electric  current  ->  Gas  breakdown  / 

Helicon  Plasma 

-  Breakdown  is  similar  to  ICP  B 

◄ — 

-  Applied  static  B0  allows  radial  and  transverse  waves 


Electric  Field 


▼  t 


▼ 


Identifying  the  helicon  mode 

-  Often  times  done  by  looking  for  axial  density  “jumps” 

•  Capacitive  (E-mode)  “pink” 

-  Hollow  Radial  Density  Profile 

•  Inductive  (H-mode)  “purple” 

-  Flat-Top  Density  Profile 

•  Helicon  (W-mode)  “blue” 

-  Sharply  Peaked  Core  Density  Profile 


More  accurately  identified  by  measuring  the  radial  density  profiles 
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Cylindrical  Helicon  Theory  (1) 


Simplified  Ohm’s  Law  neE  =  jxB0 
With  perturbations  /  ^  f0  exp  i(mQ  +  kz-  cot) 


(  f  v1 
co  juQenQ 


B 


Result  < 


yk  B0  J 


Vx5 = aB 


VxB 


where  a  = 


co  jU0en0  co  co'p 


k  Bn  k  coc A 


V2B  +  a2B  =  0 

„  1  , 

In  cylindrical  geometry  for  B_  B_  + —B.  + 

r 


Vxf  = 


dB 


dt 

VxB  =  Ju0j 

vq/  =  o 
vds  =  o 

Vx Vxg  =  V(VDg)-V2g 
VxVxg  =  -V2g 


f  2X 

2  m 


Tz 


V 


B_  =  0  where  T2  =  a2  +  E 


y 


Finite  Solution  at  r=0  Jm(Tr) 


B  =  C,J  (Tr)  B  =f3 

Z  3  m  V  X  r  rr-r  2 

1 

—aJ  -kJ, 

{r  ) 

r 

Be=  l 

U  rj-i  2 

f  /w  ,  r  r-  ^ 

—  kJm+aJm 
{r  ) 

Apply  Boundary  Condition 
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Cylindrical  Helicon  Theory  (2) 


•  Radial  Current  Density  vanishes  on  the  boundary 


j=0  @  r=a 


v  x  5  =  M(,  j 

VxB  =  aB 

MoJr  =  Bz  ikBg  =  0 
r 

im  B_  ikB0  =  aBr 
r 

Helical  Antenna 
m  =  +1 


f 


m 


aJ '  -  kJ 


V  r 


m 


m 


\ 


maJm  ( Ta )  +  kaJm  ( Ta )  =  0 


Dispersion  Relation  for  m=+l  Helicons 
Uniform  radial  density  profile 
Cylindrical  coordinates 


B0  _  eju0a  co  3.83  _  cd  n0e/u0 

n0  Zm  k  a  k  B0 


Has  been  initial  point  for  previous  work 
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Internal  Wave  Structure 


r=  5  cm 


r=  10  cm 


r=  15  cm 


B  =  100  Gauss 


B  =  250  Gauss 


B  =  500  Gauss 


B  =  1000  Gauss 
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Transition  to  spherical 


Cylindrical 

fQ  exp  i(m6  +  kz-  cot) 


Spherical 

r  i  i  x  exp  iimd)  +  kr  -  cot) 

y0  exp  i(m<p  +  kr-  cot)  vs.  j0  —  - - 

r 


Uniform  radial  density 


Assumes  no  perturbation  in  9  direction 
Waves  travel  at  a  constant  radius  V 


0  <  r  <  oo 
0  <  6  <  2n 
—oo  <  z  <  +oo 


r  +  6  +  z 


6  =  6, 


wall 


r  =  0 
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Spherical  Waves  (1) 


•  Complex  solution 

•  Requires  numerical  analysis 


B  =  B0z  =  B0  (cos  Or  -  sin  00) 
B  =  B0z  =  B0r 
oN  x  B  =  V  x  j(cos  0)  V  x  i?  j 


J  =  Jo 


exp  i{m(/)  +  kr-  cot) 
r 
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Spherical  Waves  (2) 


Separable  in  Spherical  Coordinates 


Planar  Waves  vs.  Spherical  Waves 


d_ 

dr 


r  2  dB^ 

r  - 

y  dr 


1 


a 


r2  sin#  dO 


sin# 


dB 

~d0 


d~  B 


a2  B  =  0 


r2  sin2  6  d(fr 

B(r,0,f)  =  R{r)®{ey »(*)  aVr2B  +  (ikr-l)Vx 

oo  l 

B{r,9,<f)  =  ^Yu  aiJi  (ar)Y"  i0’  0) 


( 


ikr  - 1 


A 


VxB 


lv  r  J 


=  0 


1=0  m=—l 
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Why  4  degree  angle 


•  Comparison  between  cylindrical  tubes  and  conical 

•  Investigate  wave  field  structure 

•  Density/Temperature/Potential  results  for  same  operating  parameters 


•  Small  angle  approximates  B  field  with  r  direction 

•  Tubes  in  use 
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Setup  Schematic 


•  Laboratory  Parameters 

-  Power  =  0-5  kW 

-  B-field  =  0-2  kGauss 

-  Pressure  =  0.2-20  mTorr 

-  Radius  =  1-15  cm 


-  Frequency  =  0.5-27.12  MHz 

-  Antenna  length  =  5-30  cm 

-  Gas  =  H2,He5Ax5Xe 


Electrostatic  Ion 
Energy  Analyzer 


90  GHz  Microwave 
Interferometer 


0-1200  Gauss 
DC  Magnet  Coils 


RF  Compensated 
Langmuir  Probe 


o 


Hajf-Turn 
Helicbl  Antenna 

I 


RF  EMI  Shielding 


0-100  MHz 
Waveform  Generator 


0.3-35  MHz  0-1000W 
Linear  Amplifier 
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Argon  Discharge  Photos 


*  No  Applied  B-field 


Power  <  100  Watts 

Power  ~250  Watts 

•  High  B-field  ~  lkGauss 
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Summary  and  Conclusion 


★ 


•  Dispersion  relations  which  lead  to  geometry  considerations 

•  Results  for  infinite  homogeneous  plasma  shown  lower  hybrid  frequency 
waves  with  lkGauss  B -field  in  agreement  with  a  few  experimental  data 
points. 

•  Effects  can  be  enhanced  at  low  pressure  and  through  control  of  RF  wave 
absorption 


Questions  ? 
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